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We present the current estimates of the Galactic merger rate of double-neutron-star 
(DNS) systems. Using a statistical analysis method, we calculate the probability dis- 
tribution function (PDF) of the rate estimates, which allows us to assign confidence 
intervals to the rate estimates. We calculate the Galactic DNS merger rate based on 
the three known systems B1913+16, B1534+12, and J0737-3039. The discovery of 
J0737-3039 increases the estimated DNS merger rate by a factor ~6 than what is 
previously known. The most likely values of DNS merger rate lie in the range 3 — 190 
Myr -1 depending on different pulsar models. Motivated by a strong correlation be- 
tween the peak rate estimates and a pulsar luminosity function, we calculate a global 
probability distribution as a single representation of the parameter space covered by 
different pulsar population models. We compare the global PDF with the observed 
supernova Ib/c rate, which sets an upper limit on the DNS merger rate. Finally, we 
remark on implications of new discoveries such as of J1756— 2251, the 4th DNS in 
the Galactic disk, and J1906+0746, a possible DNS system. 
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Abstract 



1 Introduction 
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Merging (or inspiralling) DNS sys- 
tems are one of the prime tar- 
gets for ground-based gravitational- 
wave (GW) interferometers such 
as GEO600, TAMA, VIRGO, and 
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LIGO (the Laser Interferometer 
Gravitational- Wave Observatory) 



( Abramovici et al. , Il992f ) . Event 



rates of the DNS inspiral searches 
by these detectors can be inferred 
using the rate estimates with an 
extrapolation out to the maxi- 
mum detection distances for any 
detector under consideration. Be- 
fore 2003, the Galactic DNS merger 
rate had been estimated between 
~ 10 -7 — 10~ 5 yr _1 (see Kalogera et 
al. 2001 and references therein). At 
that time, there were only two sys- 
tems available for e mpirical studies, 
PSRs B1913+16 dHulse fc Tavloi . 
1975[) and B1534+12 ( |Wolszczan . 
199l[ ). We calculated the PDF of 
Galactic DNS merger rates, P(R), 
based on these two systems (see 
Kim et al. 2003, hereafter KKL, 
for further discussions). Soon af- 
ter the discovery of the highly rel- 
ativistic system PSR J0737— 3039 
, we were able to revise P{R) in- 
cluding PSR J0737-3039 in col- 
laboration with the observation 



uncertainties; (iii) constraints on up- 
per limits for rate estimates based 
on the observed supernova rate in- 
corporated with our theoretical un- 
derstanding of the SN-DNS relation; 
(iv) the approximate contribution 
of J1756-2251 to the Galactic DNS 
merger rates and uncertainties in 
our assumptions on the efficiency of 
the acceleration searches for Parkes 
multibeam pulsar survey (PMPS). 
Finally, we discuss implications of 
the most recently dis covered pulsar 
binary J1906+0746 ((Lorimer et al 
2005bh . 



2 The Galactic DNS merger 
rate 



Here, we describe the main compo- 
nents of the calculation of the com- 
bined P(R) considering the three ob- 
served DNS systems in the Galactic 
disk. The merger rate of a given pop- 



team 1 ((Burgav et al.l . l2004tlKalogera et 
200"1 . The discovery of J0737-3039, 
resulted in a significant increase in 
the estimated Galactic DNS by a 
factor of ~6. This implies a boost 
in event rates for DNS searches for 
GW interferometers. In this paper, 
we present our recent results on: (i) 
the PDF of Galactic DNS merger 
rate estimates with updated obser- 
vations; (ii) the global PDF of rate 
estimates considering the systematic 
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1 Only the millisecond component 
(J0737-3039A) is considered in our cal- 
culation. For instance, the current age 
of the system is derived from the pulsar 
A. 



where A^psr,; represents the number 
of pulsars in our Galaxy with pulse 
and orbital characteristics similar 
to an observed sample i (e.g. PSR 
J0739— 3039) and /b,i is a correction 
factor to take into account pulsar 
beaming (typically ~ 6) 2 . riif e ,i is 
the lifetime of system i based on its 
observed properties. In order to cal- 
culate AfpsR,i, we distribute a large 
population of pulsars all similar to 



2 This is based on polarimetry mea- 
surements of B1913+16 and B1534+12. 
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the system i, in a model galaxy as- 
suming spatial and luminosity dis- 
tributions. Since pulsar luminosities 
are drawn from a distribution, the 
observed flux and estimated distance 
for the DNS system are not relevant 
in our calculation. Moreover, the se- 
lection effects for faint pulsars are 
implicitly taken into account. Once 
we generate a pulsar population with 
a size iVpsR, we can then calculate the 
number of pulsars detected (A^et) 
by large-scale pulsar surveys. We re- 
peat the survey simulations with a 
detailed modeling of selection effects 
for observed DNS systems. For a 
fixed iVpsR, we find that A^et follows 
a Poisson distribution, P(iV det ; iV det ), 
where N det is a mean value of iV det 
for a given population size (Apg R ). 
We require iVdet = 1, i.e., we con- 
sider only one observed system, and 
calculate the best-fit value of A" det . 
With a wide range of Apg R , we find 
-Adet = cAp SR as expected where c 
is a constant. Applying Bayes' the- 
orem to these results, we calculate 
a P(Apsr), a PDF for the popula- 
tion size of a given system i knowing 
that there is one observation. The 
calculation P(R) from P(A r PSR ) is 
straightforwad, and the full deriva- 
tion and formula can be found in 
Appendix A of Kim et al. (2004). 

The lifetime of a DNS system is de- 
fined by riife = r sd + r mrg , where r sd 
is the s pin-down age of a r e cycled 



pulsar ( Arzoumanian et al. . 1999) 



and r mrg is the remaining lifetime 
until the two neutron st ars merge 
( Peters fc Mathews! . Illffil . Based on 
the most recent observations, we es- 
timate the lifeti me of J073 7 — 303 9 
to be ~230Myr (|Lvne et all l2004h . 



This is the shortest among the 
three observed systems. The beam- 
ing correction factor f h is defined 
as the inverse of the fractional 
solid angle subtended by the pul- 
sar beam. Its calculation requires 
detailed geometrical information on 
the beam. Following Kalogera et 
al. (2001) and updating the values 
with recent observations, we cal- 
cualte f h ~5.7 for PSR B1913+16 
( Weisberg fc Tavlorl . lioolh and ~6.0 
for PSR B1534+12~ lstairs et al " 



2004) Without good knowledge of 
the geometry of J0737— 3039, we as- 
sume /b,jo737 to be the average value 
of the other two systems (~ 5.9). 

In Fig. 1, we show P{R) for a ref- 
erence model (Model 6 in KKL). 
We obtain the most likely value of 
R ~ 71 Myr -1 , larger by a factor of 
~ 5.5 than the rate estimated be- 
fore the discovery of J0737— 3039. 
The increase factor is found simi- 
lar for all pulsar population models 
we examined. The increased merger 
rates imply a boost in the inferred 
detection rate of GWs from DNS 
inspirals for ground-based GW inter- 
ferometers such as LIGO. In order 
to calculate the detection rate (D), 
we assume a homogeneous distribu- 
tion of galaxies in nearby Universe 
and a spherical symmetry in detec- 
tor sensitivity. Then, we can write 
D = i? ga i x A" ga i, where A" ga i is the 
number of galaxies in the detection 
volumne (Vd e t)- We calculate the 
number density of galaxies derived 
by the observed blue luminosity den- 
sity, (rigaj = 1.25 x lO^Myr" 1 (see 
Phinney (1991) and Kalogera et al. 
(2001) for more details). The detec- 
tion volume of LIGO can then be 
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Fig. 1. P(R) is shown on a log scale. The 
thick solid line is the Galactic rate es- 
timate overlapped with results for indi- 
vidual observed systems (dashed lines) . 
Dotted lines indicate confidence inter- 
vals for the rate estimates. The same 
results are shown on a linear scale in 
the small inset. All results are from our 
reference model. 

defined sphere for a given de- 
tection distance (20 Mpc and 350 
Mpc for the inital and advanced 
LIGO, respectively), and the num- 
ber of galaxies within Vdet is simply 
%ai x Vdct- For our reference model, 
we find that the most probable event 
rates are about 1 per 30yrs and 1 
per 2 days, for initial and advanced 
LIGO, respectively. At the 95% con- 
fidence interval, the most optimistic 
predictions for the reference model 
are 1 event per 9yrs and 1.6 events 
per day for initial and advanced 



LIGO, respectively. More details can 
be found in Kalogera et al. (2004). 

As shown in the Fig. 1, the Galac- 
tic DNS merger rate is dominated by 
PSR J0737-3039. Recently, Lorimer 
et al. (2005a) suggested that the cur- 
rent age of J0737-3039 is ~ 30 - 70 
Myr. This implies a even shorter life- 
time (riife ~ 115 — 155 Myr) knowing 
that the estimated merger timescale 
of this system is ~85 Myr. Based on 
their results, we find the most likely 
value of R for the reference model is 
~ 90 - 110 Myr" 1 . 

The beaming correction for J0737— 3039 
is not yet constrained and we assume 
MSPs discovered in DNS systems 
are not very different. As a con- 
servative lower limit, without any 
beaming corrections for all observed 
systems for a reference model, we 
obtain ~ 121" Myr" 1 with a 95% 
confidence interval. The correspond- 
ing detection rates for initial and 
advanced LIGO are 5l^ 2 x 10~3 yr -1 
and 27^21 yr -1 , respectively. Only 
when the axis geometry of J0737- 
3039 becomes available, we will be 
able to constrain the uncertainties of 
the beaming fraction, and in turn, 
the rate estimates. 



3 Global probability distribu- 
tion of the rate estimates 

In KKL, we showed that empirical 
DNS merger rates are strongly de- 
pendent on the assumed luminosity 
distribution function for pulsars, but 
not on the pulsar spatial distribution. 
Therefore, we can consider only the 
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rate dependence on the pulsar lumi- 
nosity function for simplicity. Here, 
we describe how we can incorporate 
the systematic uncertainties from 
these models and calculate, P g (R), a 
global PDF of rate estimates. Note 
that the results available on prior 
functions for pulsar luminosity dis- 
tribution are currently out of date. 
Specific quantitative results could 
change when constraints on the lu- 
minosity function are derived from 
the current pulsar sample. 

We assume a power-law luminosity 
distribution for a radio pulsar lumi- 
nosity function f(L). This function is 
defined by two parameters: the cut- 
off luminosity L min and power-index 
p. We assume prior distributions for 
these two parameters and calculate 
Pg(.ft). We fit the marginalised like- 
lihood of L mm and p pre sented by 



(|Cordes fc Chernofil Il997h and ob 
tain the following analytic formulae 
for prior functions, i.e. /(Lmm) and 
g(p): f(L min ) = a + a 1 L min + a 2 Ll lin 
and g(p) = IO^+^+^p 2 , where a { 
and Pi (i = 0, 1, 2) are coefficients we 
obtain from the least-square fits and 
the functions are defined over the 
intervals L min = [0.0, 1.7] mJy kpc 2 
and p = [1.4, 2.6]. We note that, al- 
though Cordes and Chernoff (1997) 
obtained f(L min ) over L mm — [0-3, 2] 
mJy kpc 2 centered at 1.1 mJy kpc 2 , 
we consider f(L min ) with a peak at 
~ 0.8 mJy kpc 2 considering the dis- 
coveries of faint pul sars with L 14 p n 
below 1 mJy kpc 2 (|CamiloL 12003) . 



Next we calculate P S (R) as follows: 

p s (R) 

= I P d PlL mln dL min P(R)f(L TDill )gip). 
In Fig. 2, we show the distribu- 
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Fig. 2. The global Pg(R) on a linear 
scale (lower panel) and the assumed in- 
trinsic distributions for L m in and p (up- 
per panels). Dotted lines represent the 
lower (SNl) and upper (SN\j) bounds 
on the observed SN Ib/c rate scaled 
by 5% of the observed SN Ib/c rates, 
600 - 1600 Myr" 1 (see text). 

tions of L min and p adopted (top 
panels) and the resulting global dis- 
tribution of Galactic DNS merger 
rate estimates (bottom panel). We 
find the peak value of P g {R) at 
only around 13 Myr -1 . We note that 
this is a factor ~ 5.5 smaller than 
the result from the reference model 
(R -71 Myr" 1 ). At the 95% con- 
fidence interval, we find that the 
global Galactic DNS merger rates lie 
in the range ~ 1-145 My^ 1 . These 
imply LIGO event rates in the range 
~ (0.4-60) x lO^yr" 1 (initial) and 
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~ 2 — 330 yr -1 (advanced). Since 
1997, the number of known millisec- 
ond pulsars has more than doubled, 
and therefore, constraints on L min 
and p and their PDFs based on the 
most up-to-date pulsar sample are 
urgently needed. 



4 Upper Limit of DNS Merger 
Rate Estimates. Constraints 
from Type Ib/c Supernovae 
Rates 



According to the standard binary 
evolution scenario, the progenitor of 
the second neutron star is expected to 
form during a Type Ib/c supernova. 
Therefore, the empirical estimates 
for the Type Ib/c SN rate in our 
Galaxy can be used to set upper lim- 
its on the DNS merger rate estimates. 
Based on Cappellaro et al. (1999), we 
adopt i^Nib/c ^ 1100 ± 500 Myr" 1 
considering Sbc-Sd galaxies, a Hub- 
ble constant H = 71 km/s/Mpc and 
the blue luminosity of our Galaxy 

-^B,gal = 9 X 10 9 -^B,sun- 

We calculate the fraction of SN Ib/c 
actually involved in the formation of 
DNS with a binary evolution code 



Star Track ((Belczvnski et all I2002L 

and estimate the rate ratio: 

7 = (^DNs/^SNIb/c) X 100 < 5%. 

Motivated by this result, we adopt 
the empirical -RsNib/c assuming 
7 ~ 5% and compare the value with 
the global PDF (Fig. 2). 



We note that our most optimistic 
DNS merger rate is R = 189±f§J 
Myr -1 at a 95% confidence inter- 
val (Model 15 in KKL). We obtain 



7 ~ 80% with respect to the cen- 
ter value of the empirical SN Type 
Ib/c rate (1100 Myr" 1 ) and the up- 
per limit of R at the 95% confidence 
interval (189 + 691 = 880 Myr" 1 ). 
This corresponds to 7 ~ 13% for a 
SN Type II rate, which is factor 6.1 
larger than that of SN Type Ib/c. 
In both cases, the most optimistic 
model is lower than the current em- 
pirical supernova rate estimates, but 
not really consistent with the results 
of population synthesis calculations. 
If we consider a upper limit of R at 
the 95% confidence interval from the 
global PDF, we obtain 7 ~ 13% and 
2% for the center value of SN Type 
Ib/c (1100 Myr- 1 ) and II, respec- 
tively. 



5 Implications of New Discov- 
eries to the Galactic DNS 
merger rate estimates 



Recently, Faulkner et al. (2005) dis- 
covered PSR J1756-2251, the 4th 
merging DNS in the Galactic disk 
from the Parkes Multibeam Pul- 
sar Survey (PMPS). The standard 
Fourier method failed to find this pul- 
sar and they reanalysed the PMPS 
data with an acceleration search 
(or a 'stack search' as described 
in their paper). In order to calcu- 
late the merger rate including PSR 
J1756— 2251, a detailed simulation 
is necessary to calculate the effect 
of the acceleration search with the 
PMPS. However, the approximate 
contribution of PSR J1756-2251 to 
the Galactic merger rate can be easily 
obtained. We find the total rate in- 
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Table 1 

Observed properties of tight DNS systems 



PSR name 


P s a (ms) 


P b 6 (hr) 


e c 


tL (Gyr) 


iv PSR 


B1913+16 


59.03 


7.75 


0.617 


0.37 


680 


B1534+12 


37.90 


10.10 


0.274 


2.93 


480 


J0737-3039A 


22.70 


2.45 


0.088 


0.23 


1680 


J1756-2251 


28.46 


7.67 


0.181 


2.03 


400-600^ 


J1906+0746 9 


144.14 


3.98 


0.085 


0.082^ 


300 



a pulsar spin period 

b orbital period 

c eccentricity 

d lifetime of a system 

e population size for a given DNS system 

' approximate A pop . Results depend on the survey modeling (see text). 

9 The nature of the companion is not yet clear. 9 t\^ c = r c + where Td is a 

death-time of the pulsar. 



creases by only ~ 4% due to the new- 
discovery. This is expected because 
PSR J1756-2251 can be identified 
CIS db member of the B1913+16-like 
population, which has already been 
taken into account in the calculation. 
Only future detections of pulsars 
from a significantly different popu- 
lation (compared to the known sys- 
tems), or from the most relativistic 
systems, will result in a non-trivial 
contribution to the rate estimates. 



Finally, we note the implications of 
J1906+0746 on the pulsar binary 
merger rates. This system has drawn 
attention due to the extremely young 
age of the pulsar (characteristic age 
of ~ 112 kyr; Lorimer et al. 2005b). 
If the companion is another neutron 
star, J1906+0746 would be the first 
discovery of a non-recycled compo- 
nent in a DNS system. Currently, the 
nature of the companion is totally 



unknown, and it can be either a light 
neutron star, or heavy (O-Mg-Ne) 
white dwarf. Assuming J1906+0746 
is a DNS system, we calculate its 
contribution to the Galactic merger 
rate. Because of its short lifetime 
(~82Myr), J1906+0746 can increase 
the Galactic DNS merger rate by 
about a factor 2. This implies that the 
current estimated DNS merger rate 
including J0737— 3039 can still be 
doubled! If J1906+0746 is an eccen- 
tric NS- WD system, such as J1141 
6545 ( van Kerkwiik fc K ulkarni 
| l99fl[> or B2303+46 (|Stokes et"al 
1985), it will be as important as 
J1141-6545, which is currently dom- 
inate the birthrate of eccentric NS- 
WD binaries. In Table 1, we summa- 
rize some important properties of ob- 
served merging DNS systems, which 
used in the calculation of P(R) in 
this paper. 
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6 Conclusion 



7 Acknowledgement 



We summarize our results on the 
Galactic merger rate of double neu- 
tron stars (DNS) in view of the re- 
cent discovery of PSR J0737-3039. 
We update the results based on the 
most recent observations. The global 
PDF of rate estimates incorporates 
the already known systematics due 
to the radio pulsar luminosity func- 
tion. The most likely value obtained 
from the global distribution is at 
only ~ 13 Myr" 1 , and a reanalysis of 
the current pulsar sample and radio 
luminosities is needed for a reliable 
assessment of the best fitting distri- 
bution. The DNS merger rate can be 
constrained by empirical Type Ib/c 
supernova rates, and we compare 
the supernova rate with the global 
PDF. Recent discoveries of relativis- 
ts systems, such as J1756— 2251 and 
J1906+0746 are important for the 
reduction of the uncertainties asso- 
ciated with the empirical rate esti- 
mates. As for J1756— 2251, we expect 
their contribution to the Galactic 
DNS merger rate is negligible, and 
this pulsar can be seen another sam- 
ple of B1913+16-like population. It 
is not clear whether the companion 
of J1906+0746 is a neutron star or 
a white dwarf. This system may in- 
crease the estimated merger rate of 
the relevant population by a factor 
of 2 whether it is a DNS or another 
eccentric NS-WD binary. 
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